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1.   INTRODUCTION 

In  recent  times  much  effort  has  been  spent  in  the 
search  for  narallelism  in  existing  nrop-rams.   Here,  the  con- 
cent of  data  dependence  plays  a  very  special  role.   The  ai 
or   this  thesis  is  to  study  conditions  under  which  there  i 
dependence  between  two  assignment  statements  inside  an  TV_fre 
loon  in  a  language  like  Fortran.   In  our  basic  model,  we  deal 
with  array  variables  coming  from  a  one-dimensional  array  and 
let  the  subscripts  be  arbitrary  polynomials  In  the  index  var- 
iables.  Generalization  to  multi-dimensional  arrays  noses  no 
additional  problems. 

In  section  2  we  establish  our  notation  and  e-ive  some 
useful  definitions  and  lemmas.   Our  main  result  (Theorem  3.1) 
alonr  with  some  examples,  remarks,  and  possible  extensions  is 
riven  in  section  3.   Finally,  section  *    is  devoted  to  a  detailed 
study  of  the  very  important  special  case  where  the  subscripts 
are  linear  functions. 

So  far,  no  work  has  been  done  in  the  non-linear  case, 
not  even  in  the  linear  case  with  more  than  two  index  variables. 
Cohagan  [1]  and  Towle  ([5],  ch.  5)  have  considered  the  linear 
case  when  the  number  of  variables  is  restricted  to  one  or  two. 
Towle  rives  a  very  large  set  of  necessary  conditions  in  the 
two-variable  case,  and  a  complicated  program  based  on  that  set. 
Prom  the  nractical  noint  of  view,  it  would  be  extremely  difficult 
to  extend  his  oroP-ram  to  the  linear  case  of  three  or  more  var- 
iables.  And  It  can  not  be  extended  to  the  non-linear  case  at  all. 
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The  necessary  conditions  we  have  obtained  are  valid  in  thl   general 

non-linear  case  with  an  arbitrary  number  of  index  variable:  . 
They  are  also  shown  to  be  sufficient  in  a  snecial  linear  case 
(Theorem  *J .  3 )  that  occurs  frequently  in  nractice.      Five  a 
short  and  simnle  nrorram  (based  on  those  conditions)  that  wor 
for  the  reneral  linear  case  with  any  number  of  variables,  (See 
the  Anoendix.)   And  this  oropxam  can  be  easilv  extended  to  the 
non-linear  case  with  onlv  minor  modifications.   f"ore  comments  on 
the  results  of  Ooharan  and  ,powle  are  civen  at  the  end  or  section  >\ . 


2.   PRELIMINARIES 

2.1   Definitions 

Let   r   denote  a  real  number. 

(1)  We  define   Lr  I   to  be  the  greatest  integer  not  greater  than 
r.   For  example,   [3  J  =  3,   1.3. 9J  =  3,   and   L-3.9J  =  -4. 

(2)  The  positive  nart   r   and  the  negative  part   r~"   of  r 

are  ^ivon  by 

r+  =  (!r|  +  r)/2 


and 


r"  =  (Irl  -  r)/2. 


The  basic  nronerties  of  r   and   r~   are  summarized 


.1  n  the  followinp  tv/o  lemmas. 


2 . 2   Lemma 


For  any  real  number   r, 

(1)  r+  _>  0,   r~  _>  0; 

(2)  if   r  _>  0   then   r  =  r   and   r"  =  0; 

(3)  if   r  <_   0   then   r  =  -r~   and   r+  =  0; 

(4)  r  =  r+  -  r"; 

(5)  (-r)+  =  r*,   (-r)"  =  r+; 

(6)  (r+)+  =  r+,   (r+)~  =  0,   (r~)  +  =  r",  (r-)-  =  0 


Proof.   The  proof  follows  directly  from  the  deflnltior  . 

2.3  Lemma 

Let   r,s,z   denote  real  numbers.   Tf  ri   <_   z  <_  s 


then 


+ 
-r  s  <  rz  <  r  s  . 


+ 
In  oarticular,   -r  <_   r  <_   r  . 


Proof.   If   r  >_  0,  then 

0  <_  rz  <_  rs 
and 

r  =  r+,   f"  =  0, 


so  that 


+ 
-r  s  <  rz  <  r  s. 


and 


so  that 


If   r  <_  0,   then 

rs  <  rz  <  0 


+ 

r  =  -r  ,   r   =  0, 


+ 
-r  s  <  rz  <  r  s. 


Q .  E .  D 


^hese  two  lemmas  will  be  frenuently  used  without 
always  beinn;  exnlicitly  mentioned. 


2.4  Definitions 


Let   y   denote  a  vector   (y,,  Vp,,..,  y  )   where 

y,  ,  yni....  y   are  all  integer  variables. 

1    d '  m 


(1)   For  any  positive  integer  L  we  define  a  vector  y   by 

Li 


yL 


(2)   A  product  term  (or  simply,  a  term)  in   y   1s  an  expression 
of  the  form 


ll  l2  £m 

•  1    'y2        ,ym 


where   £  ,  L.,,,.  I        are  non-negative  integer  constants. 
1    d  m 

If  for  some   r,  £   >  0   then  we  say  that  the  term  contains  the 

»   r  J  

variable   y  .   The  degree  of  the  term  is  the . non-negative  in- 

'  r  ■ 

teger      (I.    +    l„    +    ...    +    i    ).    (The    term  with   degree      0      is    the 
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constant      1.  ) 

(3)   A  polynomial  in   y   is  a  finite  sum  of  the  form 

y  ak  Pk(v")>    where  for  each   k,  aR   is  an  integer  constant 
K 

and   P.  (7)   a  nroduct  term  in   y.   The  degree  of  this  noly- 

K. 

nomial  is  the  largest  of  the  degrees  of  all  the  product  terms 
Pk(y)   such  that   ak  *   0. 

2.5  Definitions 

An  assignment  statement   S   is  a  statement  of  the  form 
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Z  =  R,   where   7.   Is  a  scalar  or  array  var^  -1 

formed  arithmatic  expression.   Z  is  cai:     e  output  variable 
of  S   and  is  denoted  by   OUT(S).   Any  variable  occurrlr 
r   is  an  input  variable  of   S,   and  the  set  of  all  input  var- 
iables of   S   is  denoted  by   IN(S). 

If   S   and   S«   are  two  assignment  statements  such 
that   S    is  executed  before  r,?i      then  we  write   Sj  <    ^?    . 

Finally,  we  have  the  followinr  convention.   Any  sum 
of  the  form  *T      (  )     where   0   1s  the  empty  set,  or  of  the 

m 
form   V  (  )      where   M  >m,  is  identically  enual  to   0. 


I 

k=M 


3.  THE  MAIN  RESULT 

We  turn  now  to  the  problem  of  data  dependence.   Con- 
sider a  program  segment  of  the  following  type: 

DO   5    I1  a  0,  n1 
DO   5     Ip  ■  0,  n2 


DO   5     I   =  0,  n 


s1(T) 


s2(T) 


sN(D 


5   CONTINUE, 
where   n>,  n„,...,  n   are  positive  integers,  T   =  (I-.,  I?'*'**    *  ^» 
and   S  ,  S  ,...,  S„   are  assignment  statements.   Take  any  two 

statements   S   and   S    (1  <  n,  q  <  N).   Let   X(f(T))   denote 
d         q      —      — 

a  variable  of   S   and   X(g(T))   a  variable  of   S  ,   where   X  is 

P  <1 

a  one-dimensional  array  and   f,  g  are  integer  valued  functions 
(such  that  the  variables  are  well  defined). 

Three  different  types  of  dependence  of   X(g(I))   on 
X(f(I))   may  arise  corresponding  to  the  following  three  possib- 
ilities : 


% 


(1)        x(f(T))  ■  OIJT(S  ),  X(pr(i))  e  II 


(2)         X(f(T))  e  IN(S  ),  X(r(D)  =  OUT<Sq); 


(3)        X(f(T))  =  OUT (3  ),  X(p,(I))  =  OUT(S„). 

For  our  purpose,  there  is  no  need  to  distinguish 

between  these  three  types.   We  simply  assume  that  one  of 

(1),  (2),  (3)   holds.   We  will  say  that  the  variable   X(r(T)) 

of  statement  S    is  dependent  on  the  variable   X(f(T))   of 

statement   S  ,   and  write 
p' 

X(f(T))  A  X(f(T)), 
if  there  exist  two  values  i      and   J  of     such  that 


(1)     Sp(i)   <   Sq(.1), 
and 

(ii)     f(i)   =  fAT). 

Our  problem  is  to  find  conditions  (necessary,  sufficient,  or 
both)  such  that 

X(f(T))  A  X(g(T)). 

In  actual  programs   f(I)   and   fr(I)   are  almost 
always  nolynomials  in   I.   We  assume  that  this  is  indeed  the 
case.   Let   d   denote  the  larger  of  the  two  dep-rees,  and 
{P,  (I)    k  £  K)        the  set  of  all  nroduct  terms  in  X  with  de- 
crees from   1   through   d.  (Here   E   is  lust  an  index  set.   If 
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the  total  number  of  these  product  terms  Is   NP,  then  we  can  take 

E  =  (1,  2,...,  NP}.)  Let 


f(I)   =   a  +  £  ak  Pk(I) 

keE 
and 


„('d    =    b  *  £  bk  Pk(I). 


kEE 

where   a,  6,  the   a's,  and  the  b's  are  all  integer  constants. 
Now  take  any  integer   L   such  that   1  <_  L  <_  m+1.   A 
polynomial  in  T  can  be  expressed  as  a  polynomial  In  TT 

Li 

(Definition  2.1\    (1))  without  any  change  in  degree.   Thus  we 


can  write 


kcE 


rL   £*_,  "Lk  lkx*L', 


where   yt   and  the  c's   are  integer  constants.   (Clearly, 

Ym+1  =  0   and   c(m+1)k  =  bk,  keE.)   Also,  let   EL1,  EL2,  EL3 

denote  subsets  of   E   such  that   ET ,   gives  all  product  terms 

in   I  (in  the  collection   {P.  (T)|k  c  E})    containing  none  of 

^t  s    IT  ,!,....  I  ;   ET  „   rrives  all  terms  containing   IT   but 

none  of   L ,,.... .  I  ;   and   ET  _   gives  all  terms  containing  at 
L+l'    '   m'        L3 

least  one  of   IT  ,,,...  ,  I  .   Then   ET , ,  ETO,  ET _   are  pair- 

L+l*    '   m  LI'   L2*   L3 

wise  disjoint  (some  of  them  may  be  empty)  and 


el1Uel2UeL3    =    e. 
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Finally,  let   n   denote  the  vector   (r^,  n?i..., 
We  have  the  following  result. 

3.1  Theorem 

If   X(f(T))  A  X(g(T)),   then 

(A)   (^  _  a)   is  divisible  by  the  gcd  of  all  the  ak ' s  and  the 
b  's,  k  e  E,  if  this  p;cd  is  not  zero; 

and  (B) 

"vl   Cak  -  cLkr  Pk(nL)  -  £  (a"  +  cLk)+  P^) 

keEL1  keEL2 

"  L   (<<  +  c^)  p,  Cn,  ) 

m<TO   k       Lk     k    l 

<_   3  -  a 

KUL1  L2 


+  y 

>-L   (a   +  c~ 


k,E   <ak  +  CLk>  W> 

''J 


for  some  positive  integer   L   such  that 


Tm+1    if   p<q 

L  <   < 

~   \  m       if   d  >_  a 


Proof    Assume  that   X(f(I))  A  X(p;(I)).   Then  there  must 
exist  two  values 
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and 


i   «   (i  ,  i2,...,  im) 


.1      vJ-ij  J  2  *  •  •  •  *  ^m 


of   I   such  that 


(1)  Sp(i)   <   S  (J),  and 

(2)  f(i)   =   K(J). 


Now  (1)   implies  that  there  exists  a  positive  integer 
L   such  that 


m  if   p  >  q; 


W  ir  =  Jr,   1  <  r  <  L  *  1  (if   L  >  1);     and 


(5) 


iL  -  ,1L  "  -1  (if   L  £  m)- 
Usinp;  the  exoressions  for   f   and   p;  we  pet  from  (2) 


that 


(6)    3  -  a  =   ^  ak  Pk(i)  "   I   bk  Pk(<1}  • 
keE  kcE 

This  implies  condition  (A)  of  the  theorem.   In  order  to  derive 
condition  (B)  we  oroceed  as  follows. 


hynotl 


I  \   VT)   "  ^,,+  I  eLk  V 


keF 
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Mence,  (6)  can  be  written  as 


(7)      3 


or   as 


- « =  -yl  + 1 1\  v*>  -  cLk  yv] 


keF 


(8)       6   -    a   =    -yL   +    £    [ak    Pk(D    -    CLk    Pk(JL)] 


keE 


LI 


+     ^     [ak     Pk(i)     "     °T,k     W1 


keF. 


L2 


+    >      ^    \(i)    "    CLk    W3 


:eE 


L3 


Now  three  cases  arise. 


Case  (i) 


1 1 '   2'** 
31,    J 2 » ' ' 


Let   k  eE   .   Then   PkU)   is  a  function  of 

.  ,  iT  ,   alone  and   p.  (TT  )    is  a  function  of 

'   L-l  k   L 

jT  n   alone.   Therefore,  (4)  imnlies  that 

lj™"  A. 


Pk(1)  "  W' 


lence,  _  _ 

ak  Pk(i)  "  CLk  W 


(ak  -  CLk}   W 
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<   (ak -  cLR)   pk(nL)        fey  Lemma  2.3, 


since   0  <  Pk(,1L)  £  pk^nL)- 


Case  (ii) 


Let   k  e  E-p.   We  can  write 


Pk(T)  .  QR(I)  Rk(IL), 


where      0    (I)      is    a   product    term   in      I      containing   none    of 

j-t  >    It.u.i.j    I    ,      and      R.  (IT  )      is   a  power   of      IT  . 
L'      L+l  m'  k      L  r  L 

Since       \(i)    =    Qk(J~L;'      We    get 


ak    VT> 
=   ak   Qk(i)    Rk(iL) 

=   ak   VV    W 

1  ak      Qk(TL)    RkUL  -   1)  by   Lemnma   2.3, 

since      0  <    VTL>   W   *  Vl5    V3L"   ^ 


<k   W 
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Hence,  a   P  (1)  -  c,^  Pk(JL) 


i  <a+k  -  cLk>  W 


i  (ak  "  CLK,+  W'   slnce  °  i  W  i  vH     ■ 


case  (iii) 


Let   k  e  E, ,.   We  have 


ak  pk(1)  "  cLk  \{h> 

iak  Pk<"L>  +  °Zk  Pk(EL>- 

since   0  <  Pk(D,   Pk(JL)  £  V"^' 

=  (ak  +  <=Zk>  W- 

The  second  part  of  condition  (B)  now  follows  fron 

(8)  if  we  use  the  inequalities  in  all  three  cases.   To  ret  the 
first  part,  multiply  (8)  by   -1   and  oroceed  as  before. 
This  completes  the  proof  of  the  theorem. 

O.F..D 

For  any  Riven  polynomial  functions   f   and   g,   we 

must  first  find  the  constants   yT  .  cT ,   and  the  sets   ET , , 

L'   Lk  LI' 

ETp,  E_  »   before  we  can  test  for  data  dependence.   Consider  the 

followinr  examoles. 


3.2  Examnle 


Consider  the  following  program: 
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DO   5    I1   -  0,  nx 


DO   5    I0  =  o,  n. 


Sd:   X(a  +  a2l?)  =  R 


S  :   X(3  +  b1il  +  b^i^)  =  r2 


5  CONTINUE. 

Suppose  that  we  want  to  test  for  dependence  of  the  output 

variable  of   S  on  the  output  variable  of   S  . 

n  p 

Here  n  =  2,  and   f(I)   and   p(T)   are  polynomials  in 

J   of  decrees  1  and   2   respectively.   We  can  take 


E  =  {1,  2,  3,  4,  c^} 


and 


P1(T)  =  llf  P2(D  =  i2, 

P3CT)     =  ij   >     VX)    =    l1l2>     VT)  =  x2 


Then   an    =   a0    =    a,,    =    a,-    =    b0    =   bn    =    br    =    0.      There   are    three 
13  4  b  2  3  5 

oossible   values    of     L:      1,    2,    3. 


Case    (i) 


Since 


Let      L   =    1.      Then      I      =    I      =    (^    -    1,    I    ) 


b1I1   +   b4IxI2 


=  b1  +  b1  d1  -  1)  +  b4i2 .+  bi((I1  -  di2, 
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we  F,et 


Yt  =  bi  »   cn  =  bi»  cio    =    hn> 


1     1»    11 


12     H 


C13  =  °»  Cl>\    =    bH»   Cl^  =  °* 


Also,   Eia  =  0,   E12  =  (1,  3),  E   -  {2,  H,  5). 


case  (11) 


Let   L  = 


2.   Then   I.  =  I2 


=  (Ilf  I2  -  1).   Since 


b1T1  +  b,!^ 


■  <bi  +  V1!  +  Vl(I2  "  1)' 


we  eet 


Y2  =  0,   c21  =  bx  +  bi|}   c22  =  0, 


c23  ~   '   °24  ~   4'   C2S 


Als 


o,   E21  =  (1,  3},  E22  =  (2,  4,  5>,  E23  =  0, 


Case  (ill) 


Let   L  =  3.   Then   IL  =  I  =  (I    I  ).   V,Te  have 


Y3  =  °'   C31  =  bl*   °  32  =  °' 


c33  ~  °'   °3^  ~  b^a   °35  ~  ° 


Also,   E31  =  {1,  2,  3,  4,  5>,  E32  =  0  ,  E33  =  0 
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Corresponding  to  the  three  values   1,  2,  3   of  L, 
the  inequalities  in  condition  (D)  of  Theorem  3.1  take  the 
following  three  forms: 

(1)    (L  =  1) 


b~   -  b,  n,  -  a"  n^  -  b^  n-\n2 


<  B  -  a 


<_  -  b    +  b~  n   +  a«  rip  +  b^  run,,, 


(2)    (L  =  2) 


-(b1  +  b^)  n1  -  a"  (n   -  1)  -  b2|  n1(n2  -  1) 


<  3  -  a 


<_  (b-j^  +  bi4)~n1  +  a?  (n2  -  1)  +  b^  n.(n2  -  1), 


(3)    (L  =  3) 


-b.  n,  -  a?  n2  -  b^  ninp 


<  6  -  a 


<  b,  n,  +  a,,  n~  +  b7  ri.n~ 
—      11     2.      2.  4   12 


Now  sunoose  that  the  variable   X(3  +  b  I,  +  b^I-Ip) 

of  statement   S    is  dependent  on  the  variable   X(a  +  a0I0) 
q        f  2  2 

of  statement   S  .   Then,  (3  -  a)   must  be  divisible  by  the  n;cd 


oT   a~,  b  ,  bi(.   Tf   p  <  q,   then  one  of   (1),   '   ,  '  • 
hold.   Tf  n  >  o,   then  one  of   (1)   and  (2)   mm 


3.3  Example 


Consider  the  followinr  nrofram  serment: 


DO   5     T1  =  0,  n1 


DO   5     I2  -  0,  n2 


Sx:   Z  =  X(p,  -  I*  +  5lJ  +  2I1I2^ 


S2:   X(«  +  I1  -  3I2)  =  R 


5  CONTINUE 

Suppose  that  we  want  to  test  for  dependence  of  the  inout  var- 
iable of   S.   on  the  output  variable  of   S  . 
1  2 

Here, 

m=2,   p=2,   q  =  1 , 

f(I)  =  a  +  I1   -  3I2, 

and  «<T)  =  3  -  I*  +  5I?  +  21^  . 

We  can  take 

E  =  {1,  2,...,  91, 

and 
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P1(I)    =    I1,    P2(I)    =    I2, 


P3(I)    =    I*,    P2](D    =    1^2,    P5(D    =    l\ 


-   T3 


T^      =      T3 


P6(I)      =      Ij>,      P7(I)      =     IJI2,      Pg(l)      =     I1I|,      Pg(l)      =     I|. 


Then, 


a,    -    1,      a g   -   -3, 


a3        ai4        a5        ^f,   "   a7   ~   a3   ~   °s 


a,,   =   a0   =   a.,    =   ar.   =   a,-   = 


b3   =   ~lj      b6   =    5»      b8   =   2' 


bl   =  b2  =  b4   =  b5   =   b7   =   b9   =   °' 


Since   p  >  q,   there  are  only  two  possible  values  of   L:   1,  2 


Case  (i) 


Let   L  =  1.   Then   IL  -  I.  =  (I-L  -  1,  I2). 


Since 


-  I?   +    51?  + 


!I,I 


1    "~1         --1 

=  -  [(I-l  -  1)  +  l]2  +  5[(I1  -  1)  +  I]3  +  2[(I1  -  1)  +  1]I2 

=  i\   +   13(I1  -  1)  +  14(I1  -  l)2  +  2I2  +  5(IX  -  I)3 


+  2(1..  -  1)1* 


we  n;et 


Tn  =  ^>   On  =  13,   c,Q  =  14, 


"11 
Jl6 


=  2,   c,A  =  5,   clR  =  2, 


13 

18 


c12    cl4 


c17    c19 


0. 


Also,   K,,  =  0,   I 


!12  =  (1,  3,  6),  E13  =  {2,    fc,  5,  y,  B, 


Case  (ii) 


since 


Let   L  =  2.   Then   T.  =  T.  =  (I_,  I-  -1). 

L        C  id 


2  1  2 


=  "  Jl    +    5I1  +  2I1[(I2  "  1}  +  1] 


.2 
*1 


=  21,  -  i,  +  'a^ip  -  D  +  51J  +  2i1d2 


i)2, 


we  ret 


Yp     >  ^21        ^2^    — 
coh=^»   Go^=5,   c0o    2, 


2^ 


26 


28 


c22   c25   C27    °29    °* 


A 


lso,   E   =  (1,  3,  6),  E22  =  {2,  A,  5,  7,  8,  9>,  E23  =  0. 


Corresponding  to  the  two  values  of   L,  the  inequalities 
in  condition  (R)  of  Theorem  3.1  take  the  followinp  two  forms: 


(1) 


(L  =  1) 


(2) 


2    c  3 
nl  "  5nl 


2 
2n,n 


-  3n^  <  6  -  a  <  -4, 


(L  =  2) 


3  +  n1  -  9n?   +  9n2  -  5n?  -  2n_Ln2  -  3n^  < 


-  a  <_  n,, 


If  the  input  variable  of  S.  is  dependent  on  the  out- 
put variable  of  S?,  then  by  Theorem  3.1  either  (1)  or  (2)  must 
hold. 


21 


3.^1  Remarks 


(1)   The  conditions  of  Theorem  3.1  are  not  sufficient 
for  data  dependence  in  general.   We  show  this  by  an  example. 

Let 

m  =  1,   n1  =  2, 

f(I)  =  2I1,   g(T)-  H  -  3IX. 


Since   d 
Then 


1   we  can  take   E  =  {1}   and   P,(I) 


a    =    0, 


ax  =  2,   (5-4,   b1  = 


=  T. 


-  3. 


Let   n  >  q.   Then   L  =  1.   Clearly, 


y1   =  -3,  cn  =  -3, 

E1n    0,   E10    (!},   E13  =  0. 


11 


12 


Condition  (R)  of  Theorem  3.1  becomes 


x 

or, 


(al  +  Cll)+(nl  "  1)    -   3  ~  a  -  "Yl  +  (al  " 


Cll^  ^nl~  "^ 


3-0  <  ^  <  3  +  5 


which  is  true.   Condition  (A)  is  also  satisfied.   However,  two 
values   i,  ,  ,j    of   I    can  not  be  found  such  that 


and 


0  <  11  <  ^  <  2 


21.,  =  H   -    33  ± 
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Therefore,  the  variable  '/.( <\   -   31.) 

dependent  on  the  variable  /.(?!) 

(2)  Suppose  for  a  moment  that   1  ,...,  i    and   .},,...,  ' 
denote  real  variables,  and  define  a  real  valued  function 

P(ilf....  im,  .1^...,  Jm)  =  V  Eak  Py(T)   -  bk  %(!)]. 

keE 

For  a  riven  value  of   L,   v/e  consider  a  subset       -  P  the 

Li 

2n-dimensional  Kuclidean  space,  defined  by 

0  <  i  ,  ,1   <  n    (  1  <  r  <  m) , 

and  the  conditions  ( ;0 ,  (5)  in  the  proof  o^  Theorem  3.1. 
Since   D,   is  connact  and  connected  and   F   is  continuous, 
the  image  set   F(DT  )   is  a  finite  closed  interval.   (T,ee 
Dieudonne'  [2],  mheorems  3.17.10  and  3.19.8.)   Condition  (B) 
of  Theorem  3.1  is  necessary  and  sufficient  for   (ft  -  a.)      to 
lie  in  the  interval   F(DT),  i.e.  for  the  existence  of  a  "real" 
noint  in   DT  at  which  the  value  of   F   is   (B  -  a). 

(3)  Suppose  that  condition  (B)  of  Theorem  3.1  is  satisfied 
for  at  least  one  value  of   L.   Let   Ln   be  the  largest  value 
of   L   for  which  it  is  satisfied,  and  suppose  that   L~  <_  m. 
Then,  there  will  be  no  data  dependence  if  the  index  variables 
I,,  T0,...,  I    are  kept  fixed  at  any  set  of  permissible 
values . 

(l\)      The  whole  idea  behind  Theorem  3.1  is  to  ret  the  sharpest 


23 

possible  bounds  for   (p  -  a)   that  can  be  directly  obtained  using 

the  nroperties  of  the  different  terms  in   f(T)   and   g(T). 
Let  us  comnare  the  second  inequality  in  condition  (B)  against 
two  other  inequalities  that  can  be  derived  much  more  easily. 
From  (6)  in  the  Droof  of  Theorem  3.1,  we  get 


-  «  <  ]T  Ca*  +  bk)  Pk(n). 


keK 


Let      m   =    1,      f(I)    ■   a   +   I,,      and      g(T)    =    8   +  -I-.      Then,    this 


will   i^ive      3  -  a   <  n.  . 


But  the  second  inequality  in  condition  (B)  gives 


0  -  a   <_  -  1     for   L  =  1 

and 

B  -  a   <  0      for   L  =  2. 


From  equation  (7)  in  the  nroof,  we  get 

8  "  n     i  "YL  +I  (ak  +  cZk>  Pk(nL). 
keE 

The  second  inequality  in  condition  (B)  is  clearly  better  than 
this.   Because,  for  any  two  real  numbers   a   and   c   we  have 


(a  -  c)   <  (a  -c)   <a  +  c  , 


and  these  relations  may  be  strict. 

One  could  give  a  similar  discussion  about  the  first 
ineouality  in  condtlon  (B). 
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'U[j    Kx tons  Ions 

(1)   Let   I   denote  any  one  of  the  indie    I.,  I,, 
...  I  .   For  simplicity,  we  have  assumed  that  I  runs  thro, 
an  ordered  set  of  values  of  the  type   (0,  1,...,  n},       n 
is  a  positive  integer .   However,  one  frequently  has  a  statement 
like 

"DO   5   I   =   M0,  nx,  e" 

so  that  I  runs  throurh  the  values   {N.,  M_  +  e,...,  .'.'}. 
(Here,  H~,  H  ,   and   e(/  0)   are  integers  and   e   divides 
I!   -  M- . )   In  this  case  we  can  introduce  a  new  index  variable 
If   by 

Tf   =   (I  -  !I0)/e. 

Then  the  above   DO  statement  mav  be  replaced  bv 


"DO   5    i'  =  0,  n", 


Where   n  =  (N,  -  U~)/e  if  every  occurrence  of   I   inside 

i 

the  loop  is  replaced  by   Ii   +  el  . 


(2)   So  far,  we  have  worked  only  with  one-dimensional  arrays. 
Suppose  now  that   X   is  multi-dimensional.   Let   ND   denote  the 
number  of   dimensions  of   X,   and   ^B^.»  LB,   the  upper  and 

lower  bounds,  respectively,  of  the   k    dimension  of   X.   Then 

i 
X   can  be  mapped  onto  a  one-dimensional  array   X    (in  a  one 

to  one  fashion)  by 


X(£  -,...,    H,m  )   *>   X    (P(£-,...,    *,»Tr\  )  ) 
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where        a    is    defined   below: 

aU,,...,    *,JD)    =    £\l\   ~   LBk) 

k=l 
with 

'(UBk+1  -  LBk+1+  1)  ...  (UBND  -  LBND  +  1)   if   k .<  ND  -  1 

if     K  =  ND. 

Thus  Theoren  3.1  can  be  applied  to  multi  dimensional  arrays 
after  they  have  been  transformed  this  way. 

(3)   Theorem  3.1   can  be  extended  so  that  the  statements 

S    and   S    are  not  required  to  be  within  the  same  loon. 

n         q 

Let  us  illustrate  this  bv  the  following  examnle: 


DO     5      I1  =  0,  n. 


DO     5      It  =  0,  n1 


S 


CONTINUE 


DO     7      It+1  =  0,  nt+ 


Sq 


DO    7      I   =  0.  n 


7     CONTINUE 
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Here,   f   Is  a  function  of   I,,...,  I    alone  and   r   '-a 

function  of   I^.t,....  I    alone.   Equation  (6)  In  the  proof 
t  +  1 '      m 

Theorem  3.1  holds,  where   1   and   j   are  now  any  two  arl  Ltn 
values  of   I.   We  ret  the  condition 

-X  (<\ +  bk>  pk("'  i  b  -  °  -  T  ut. +  bk>  pk(n_) 

keE  keE 

which  replaces  condition  (R). 
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4.   THE  LINEAR  CASE 


In  this  section  we  assume  that  both   f(T)   and   p5(T) 
are  polynomials  in   I   of  degree  1.   This  is  the  most  important 
special  case  that  one  encounters  in  practice.   Our  first 
theorem  here  is  the  linear  version  of  Theorem  3.1. 


4.1  Theorem 


Let 


m 


f(T)  =  «  +  y  akIk 


and 


m 


„(i)  -  u  ♦  y  bkik 


where   «,  3,   the  a's,   and  the  bfs   are  all  integer  constants 


If   X(f(I))  A  X(g(T)),   then 
(A)   (8  -  a)   is  divisible  by  the  ";cd  of  all  the  a.  '  s 

and  the  b.  's,   1  <  k  <  m   if  this  o;cd  is  not  zero: 

k   '     —   —  >  ' 

and      (B) 
L-l 


m 


-bL  .  y   (ak  -  bk)  nR  -  (a"  +  bL)  (nL  -  1)- 


(a!~  +  b,  )n, 
k    k  k 


y 


l 


k=?L+l 


<  8  -  a 


L-l 


-  -bL  + 


fak  "  bk)+nk  +  (alt  "  bL)  +  (nL 


-  1 


)  +  £  (a*  ♦  bj)n, 


=  1 


k=L+l 


for  some  oositive  integer   L   such  that 


L     < 


m    +    1  ir  D    <    fl 


2fl 


n 


if         p    >    q  . 


(For      L    =    rn    +    1,      we    define      a,     =    b.     =    n,    =    0.) 


Proof  We   will    use    the    notation    of    section    3.       Since      d  ■   1, 


We   can   take 


and 


Case    (i) 


'hen 


=    { 1 ,    2  ,  .  .  .  ,    rn } 


Pk(I)    =    Tk 


,       1    <    k    <    m, 


Let         1    <    L    <    m. 


Ej.    -       {1,    2,...,    L   -    1}, 


F,L2    =       (L), 


ET  _    =       {L   +    1,    L   +    2,  .  .  .  ,    m}. 


(If      L   =    1      then      ET  ,    *    0f.      IP      L   =   m      then      ET  -    =    #.  ) 


We   have 


I 

keF 
m 

I 

L-l 


bk    Pk(I) 


b,  I, 
k   k 


m 


Zb,  I,  .  +   bT  IT     +      V 
k  k|      l  l       ; 

K=l  k=Wl 


b.  I. 
k   k 


-  \ +  £ 


m 


k 


k=l 


b,  I.     +   bT  (Ir    -    1)    +     V    b.  L 

k    k  L      L  /         k    k 

k=l+L 
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=  bL +  y  \  W 

keE 


'herefore, 


and 


\  '  \ 


cT1     =   b,  ,    k   e   E 


'Lk  k 


.*> 


case    (ii) 


mhen 


Let      L  «  m   +    1. 


EL1=    E      and     EL2    =  EL3    =  0. 


.Since   IT  =  I,   we  <ret 


L 


Ibk  Pk(I)  =  Ibk  W' 
keE  keE 


so  that 


and 


yl  ■  °  "  \ 


c  T .     =   b ,        .  k  e  E 

Lk     k        ' 


rphe  conditions  (A)   and   (B)  of  the  present  theorem 
now  follow  from  the  similar  conditions  of  Theorem  3.1. 

Q .  E  .  D . 

Theorem  4.1  is  easier  to  anply  than  Theorem  3.1  as 
we  see  below. 


V 


H.2    Kxamole 


Let 

m  •  3, 


and 


Then 


f(I)  -  o  +  HI1  -  Bl2, 
r(T)  =  B  -  6I1  +  4l2  +  161 


ai  =  4»   a2  =  ~8,   a3  =  °» 
b -  -  -6,   b2  =  ^,   b   =  16. 

The  red  of  the  a's  and  the  b's  is  2.   Corresoondinp  to  the 
possible  values   1,  2,  3,  ^   of   L,   the  inequalities  in  con- 
dition  (B)   of  Theorem  *1 . 1  take  the  following  rorrns : 

(1)       (L  =  1) 


6  -  12n2  -  l6n    <  3  -  a  <_  -4  +  10n1, 


(2)  (L  =  2) 

8  -  12n2  -  l6n   <  B  -  a  <_  -*i  +  lOn  , 

(3)  (L  =  3) 

-12  n2  -  l6n-  <_  B  -  a  <_  -16  +   lOn..  , 

CO       (L  =  1|) 


-12n2  -  l6n„  1  B  -  a  <_  lOn 
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Suppose  now  that    X(f(I))  A  X(r,(T)).   Then  by 

Theorem  '1.1,  (8  -  a)   must  be  divisible  by  2.   If  p  <  q,   then 
at  least  one  of   (1)  -  (4)   must  hold.   If   p  _>  q,   then  at  least 
one  of  (1)  -  (3)   must  hold. 

A  computer  program  that  tests  for  data  dependence 
usinr  the  conditions  of  Theorem  4.1  is  riven  in  the  Appendix. 

The  example  in  Remark  3.4  (1)   shows  that  the  conditions 
of  Theorem  4.1  are  not  sufficient  in  the  peneral  linear  case. 
Fortunately,  they  are  sufficient  in  a  special  linear  case  that 
occurs  frequently  in  actual  nroprams.   This  is  shown  by  Theorem  4.3 


4.3   Theorem 


Let 


m 


f(I)  =  a+  t£   ukIk 
k=l 
m 
and  r(D  =  3  +  tV  vkIk, 

k=l 

where   ct,  8,  t   are  integer  constants,   t  >  0,    and 
u,  ,  v,  e  {-1,  0,  1  •  },  1  _<  k  _<  m. 
Then,   X(f(D)  4X(g(!))   if  and  only  if 

(A)   t   divides  ( 8  -  a) , 

and  (B) 

L-l  m 


Y    (uk  "  vkr\  "  {ul   +  vr/(nL  -  1}  -    ^    (uk  +  vk)ni 

k=l  k=L+l 


u 


<    (3  -  a)/t 


L-l 


m 


<  -v, 


+  I   (Ul 


"  V\  +  K  "  vl/(nL-  1}  + 


+  7" 


k=l 


I 


for  some  positive  integer   L   such  that 


L  < 


m+1    if    n  <  n 


m 


if   p  >  a 


(For   L  =  n  +  1,   We  define   uT  =  vT  =  n,  =  0.) 
'  L     L     Jj 


Proof    Since  the  'only  if  nart  follows  directly  fron  ^heorem 
4.1,  we  need  prove  only  the  'if  part.   Assume  that  the  conditions 

(A)  and  (13)  hold. 

Let  (0  -  cx)/t  +  v,  >  0.   (The  negative  case  can  be 
handled  in  an  exactly  similar  way.)   The  second  inequality  in 

(B)  imnlies  that 


0   <_    (3   -    ot)/t    +    V! 

L-l 


m 


1      £      (uk-   V+nk    +    (UL    -    VL)  +  (nL    "    1}    +    X         ^    +   ^)rV 


k=l 


k  =  L+l 


Therefore,    we    can   write 


m 


(1)       (P   -    a)/t    +   vL   = 


I  v 


k=l 


where 


if 


(2)      0    <    hR    <_ 


(Kk    "   vk)+nk 

(u      -   vT)    (nT     _    i) 


Cu*   +   Vk)nk 


1    <    k    <_  L   -    1 
L    -    -w         if      k   =   L 

if      L  +   1   <  k  <  m. 
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Now,  in  order  to  show  that   X(f(I))  A  X(g(I)), 

it  will  suffice  if  we  can  find  two  values  i   =  (in i  ) 

1'    ■   rn 

and  J   =  (J  , .  . . ,  ,1  )   of  T   such  that 


(3)    S  (i)    <   S  (j) 
v       p       o    q 


and 


('0   f(i) 


g(T). 


Substituting  for   f   and   g   in  (4)  we  get 


m 


m 


a  +  t   r  ukik  =  B  +  t  y  vkJk 


m 


or 


Z[- 


M 


i   _  v,  ,1,  ]   =  (B  -  a)/t. 
k  k    k'  k 


k^ 


k=l 


By  virtue  of  (1),  this  is  equivalent  to 


m 


(5) 


I   [Vk  "  Vk]  +  ^L1!  "  VL('1L  "  1)]=  I 


V 


k=l 
k^L 


k=l 


We  will  find   i   and  J  such  that  (5)  and  (3)  are 
both  satisfied.   This  will  be  done  by  picking  suitable  solutions 
to  the  individual  equations 


(6) 


Vk  -  Vk  =  hk;  1  -  k  -  m'  k  *  L; 


uTi-  v,  (j   -  1)  =  hT   (if   L  <  m). 


We  give  below  a  table  that  lists  these  solutions.   For  any  given 

k,   detailed  values  of   u,   and  v   are  considered  only  when 

k        r 

h.  j*  0.   If  h,  =  0,   then  we  always  take  the  trivial  solution. 
Remember  that  the  range  of   h    is  given  by  (2). 


Range  of  k 

Uk 

vk 



Ranpe  of  h 
k 

Equation 

Solution 

he 

•>k 

l<k<L-l 

hk-o 

Vk-Vk"0 

0 

0 

0 

-1 

0<hk<nk 

^k  ■  hk 

0 

hk 

(if  L>1) 

1 

0 

0<hk<nk 

lk  ■  hk 

hk 

n 
k 

1 

-1 

0<hk<2nk 

Whk 

Lhk/2j 

hk-Lhk/2j 

k-L 

hL  =  0 

Vl-vl(Jl"1h 

0 

1 

-1 

-1 

°£hL£nL-l 

"W^L 

0 

V1 

(if  L<m) 

0 

-1 

0<h. <nT -1 
—  u—  l 

h-1   '   hL 

0 

V1 

1 

0 

0£hL<nL-l 

hmhL 

hL 

nL 

1 

-1 

0<hL£2n.-2 

W^L 

LhL/2j 

hL+l-LhL/2j 

L+l<k<m 

v° 

"kWk"0 

0 

0 

1 

1 

0£hk<nk 

1k--'k-hk 

hk 

0 

(if  L<m) 

-1 

-1 

0£hk<nk 

-V^k-^t 

0 

hk 

0 

-1 

0<hk<nk 

Jk-hk 

0 

hk 

1 

0 

0£hk<nk 

Vhk 

hk 

0 

1 

-1 

0<hk<2nk 

Vrt 

Lhk/2j 

hk-Lhk/2j 

Thus  two  values  T  =  (i, ,...,  i  )  and  J  ■  (J-,,...,  J„) 

l       m  l       m 

of  T  can  be  found  such  that  (6)  is  satisfied  and 


0  <  \   <  Jk  <  \»        1  <  k  <  L  -  1   (if  L  >  1); 
0  <  iL  <  JL  -  1  <  nL  -  1,         (if  L  <  m); 

0  <  1k»  Jk  -  nk»  L  +  x  I  k  i  m       (if  L  K   m)« 
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This  means  that   i   and  j   satisfy  (3)  and  (4).   Therefore, 

X(f(T))  A  X(g(T)). 

Q.E.D 

Let  us  illustrate  the  application  of  Theorem  4.3  by 
an  example. 

4.4  Example 

Consider  the  following;  program  segment: 

DO   5    I1   =     0,  10 
DO   5    I?  =  0,  5 

DO   5    I3  =  0,  9 

Sx:   X(13  +  6l1  -  61  )  =  R1 


Spt   X(l  +  61   -  6l2  +  61  ')  =  R2 


5      CONTINUE. 

And  suppose  that  we  want  to  test  for  dependence  of  the  output 
variable  of   S?   on  the  output  variable  of   S, . 

Here  we  have 

m  =  3,   ni  =  10>  n2  =  5,  n„  =  9, 

p  =  1,  q  =  2, 

f(T)  =  13  +  6I1  -  6I3 


g(I)  =  1  +  6l1  -  6lp  +  61-,. 
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f(T)   and   g(T)   have  the  forms  needed  for 

a  -  13,  p  -  1,   t  =  6, 

u1  =1,   u2  =  0,   u^  =  -1, 

vx  =1,   v2  =  -1,  v^  =  1. 

t   divides   (fl  -  a).   And  for   L  =  1,  the  inequalities  in  condition 
(13)  of  Theorem  Jl .  3  become 


"Vl  "  ^Ul  +  Vl^  ^nl 


1)  -  (u~  +  v2)n2  -  (u~  +  v^)n^ 


<  (3  -  «)/t 


<_  -  v1  +  (u1    -  v^)    (n1 


1)  +  (u2  +  v~)n2  +  (u_  +  v^)ru 


or,     -28  <   -2    <    k 


which  is  true.   Therefore,  by  Theorem  H . 3  the  outnut  variable 
X(l  +  61   -  6lp  +  61-,)   of  statement   S?   is  denendent  on 
the  outnut  variable   X(13  +  61,  -  6l_)   of  statement   S, . 
Indeed,  taking 


i  =  (0,  0,  0) 


and 


.1  =  (1,  0,  1)  , 


we  see  that 


S,  (i)  <  S0(.1) 
1     o   2 


and 


f(D  =  g(J). 
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4 . 5  Remarks 

We  g-ive  now  some  comments  on  the  results  of  Cohap;an 
[1]  and  Towle  ([51,  ch.  5). 

Theorem  1  of  Cohafran  can  be  derived  from  Theorem  4.1. 
His  Theorem.  2  and  Lemmas  2.1,  2.2  follow  immediately  from  well 
known  results  in  the  theory  of  Diophantine  Equations.   (See 
Mordell  [4],  ch.  5,  Sec.  1.) 

Now  consider  Towle.  Parts  (1)  and  (3)  of  his  Theorems 
5.3  and  5.4  can  be  derived  from  our  Theorem  4.1.   However, 
part  (2)  of  each  of  these  two  theorems  is  false.   (Take  n,  =  10, 
p  =  20,  q  =  4,  a,  =  1,  b2  =  2.   Since   20  +  0  =  4  +  2*8, 
there  is  data  dependence.   But   (d  -  q)/(b   -  a,  )  =  16  >  n  ) . 
His  theorems  5.5  and  5.6  are  contained  in  the  nroof  of  our 
Theorem  3.1.   Towle  uses  a  theorem  of  I7an([3],  Th .  10)  on 
the  inconsistency  of  a  system  of  linear  equations  to  derive 
a  set  of  necessary  conditions  for  data  dependence.   These  con- 
ditions should  be  equivalent  to  the  inequalities  in  (B)  of 
Theorem  4.1,  since  both  sets  are  (separately)  necessary  and 
sufficient  for  the  existence  of  real  solutions.   However, 
these  conditions  are  rather  imnractical.   Because,  they  are 
difficult  to  derive,  and  even  in  the  two  variable  case  with 
L  =  1,   there  are  29  of  these  conditions  for  the  sinple  pair 
of  inequalities  in  (B)  of  Theorem  4.1.   (See  Towle  [5], 
Appendix  B. ) 
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APPENDIX 

Here  we  reive  a   PL/I   program  that  tests  for  data 
dependence  usinp:  the  conditions  of  Theorem  ^.1. 

The  innut  to  the  program  consists  of  the  values  of 
W,  ALPHA,  A,  BETA,  B,  N,  P,  Q   where 

M  =  m; 

ALPHA  =  a ; 

A   is  an  array  of  size   M  with   A(K)  =  a,, 

1  <  k  <  M; 
BETA  =  3; 
B   is  an  array  of  size   M  +  1   with   B(K)  =  b.  , 

1  <  k  <  M,  B(M  +  1)  =  0; 
N   is  an  array  of   size   M  with   M(K)  =  n,  , 

1  <  k  <  M; 

P  -  p,  Q  -   q. 

The  gcd  of  the  set   (A(l),...,  A(M),  B(l),...,  B(M)} 
is  comnuted  and  condition  (A)  of  Theorem  '4.1  is  tested.   If 

it  holds  then  the  program  goes  on  to  test  condition  (B).  ((A) 

i         i  it 

is  referred  to  as  the   GCD  Test  ,  and  (B)  as  the   Minimax  Test  .) 

Mote  that  condition  (B)  can  be  written  as 
-PSI(M,L,  -A,  -B,  N)  <  BETA  -  ALPHA  +  B(L)<  PSI(M,  L,  A,  B,  M), 

where 

L-l 
PSI(M,  L,  A,  B,  M)  +  y-j  (A(K)  -  B(K))  N(K) 

+  ((A(L))   -  B(L))+(M(L)  -  1)  +  £  ((A(K))+  +  (-B(K) )+ )N(K) . 

k=L+l 


Since  PS]  Ls  always  non  negative 
of  the  two  inequalities  in  the  condition,  d< 
of  [BETA  -  ALPHA  +  B(L)].   The  number  of  tines  the  fimctJ  I 
Is  comoutod   (and  then  used  in  a  comparison)  ls  at  mo^t   N  ♦  1, 
The  similarity  between  the  conditions  o?   Theor'     .1 
and  ^  .  1  shows  that  this  prop-ram  can  be  extended  to  the  non  lin' 
case  without  too  much  difficulty.   The  only  orobl<-      that  we 
must  know  the  sets  BL , ,  ^Tp>  Ej  ^   and  the  interers  y.  ,  c. . 
and  incornorate  then  into  the  program. 
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PSI  ;  PROCEDURE  (M.L.A,B,N)  RETURNS  (FIXED  BIN  (31)): 

PLUS  :  PROCEDURE  (R)  RETURNS  (FIXED  BIN  (31)); 
DCL  R  FIXED  EIN  (31)  ; 

IF  R>0  THEN  RETURN  (R)  ; 

It  SB  F.ETUKN  (0)  ;       


END  PLUS^; 

DCL_XM.  L.  A(*),  B(*).  N(»),  K,  SUM)  FIXED  BIN(31); 
3TJH  =  ~^~ 


DO  K  =  \    TO  L-1 ; 

SUM  =  SUM  +  PLUS  (A  (K)  -  B  (K)  )  *N  (K)  ; 

rl^L<=M  THEN  SUM  =  SUM  ♦  PLUS  (PLUS  (A  (L) )  -  B  (L)  )  *  (N  (L)- 1)  ; 

DO    K  =    L+1    TO   H; 

SUM   =    SUM    ♦     (PLUS  (A  (K)  )     ♦    PLUS  (-B  (K)  )  )  *N  (K)  ; 

END: „ 

PUTSKIP   LI3T~f •  PSI=',  SUM)  ; 

RETURN     (SUM) ; 
END   PSI; 

HCL    (M,    Lr    ALPHA,     BETA,    P,    Q)     FIXED    BIN131)  ; 

DCL     \K.    DIF,    GCD,     X,     3,    UB)     FIXED    BIN(31); 
DCL    DEPENDENT    BIT  (1)  ; 
_9JL  ENDFILE  (SYSIN)     STOP; 

DO   V'HILE     ('  1'B)  ; 
GET    LIST     (M)  ; 

BEGIN; 

DCL     (A(M),    B(M*1),    N(H)J     FIXED   BIN(JI); 

GET    LIST     (ALPHA,     A,    BETA-    B,    N,    P,    Q)  ; 
PUT    SKIP(3)     DATA     (M,    ALPHA,    A,    BETA,     B,    N,    P,    Q)  ; 

GCD   =    AM);  /*    COMPUTE    THE   GCD    OF    ELEMENTS    OF    A    S    B   */ 

DO    K    =    1    TO    M-1 ; 
X   =    AJK+1) ; 
DO   WHILE     (X    -.=    0)  ; 

R    =    MOD(GCD,    X)  ; 
GCD   =    X; 


X    =    F; 

EWj 

END; 

DO    K   =    1    TO    M; 


DO   WHILE     (X   -i=    0)  ; 

R    =   MOD    (GCD,    X)  ; 
GCD   =    X; 
X   =   R: 


END; 
END; 
PUT    SKIP    DATA     (GCD) ;  /*    GCD    HAS    BEEN    COMPUTED    */ 


D I F   =    BETA    -    ALP  H  A  ;' 

IF    GCD   -.=    0    S    MOD  (DIF,    GCD)     -»=    0    THEN 

PUT    SKIP    LIST     (•NOT    DEPENDENT    BY    THE    GCD    TEST1); 
ELSE    DO: 

IF    P    <    C    THEN    UB    =   M+1;    ELSE    OB   =   ft; 


DEPENDENT    =    •  O'B; 
DO    L    =    1    TO    UB    WHILE     (-DEPENDENT) ; 
K    =    DIF    ♦    B(L)j 


I'F     (K>0"6"K    <=    PSI  (M,L,A,B, N) ) 

|  (K<0    S    K    >=    -PSI  (M,L,-A,-3,N)  )  |  (K=0)     THEN    CO; 
DEPENDENT   =    • 1 • B; 
PUT    SKIP    LI ST     ('PROBABLY    DE PENDENT' .' 1=' . L) ; 


I  END; 

IIP   DDT; 


END; 
END  * 
IF   -DEPENDENT    THEN    PUT    SKIP    LIST 

('NOT    DEPENDENT    BY    THE    MNIMAX    TEST'J_; 

"END;  ■" 

END; 
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